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Abstract 

By means of Monte Carlo simulation with local spin update Metropolis algorithm, we have elucidated non-equilibrium 
phase transition properties and stationary-state treatment of a disordered binary ferromagnetic alloy of the type ApBi^p 
on a square lattice. After a detailed analysis, we have found that the system shows many interesting and unusual 
thermal and magnetic behaviors, for instance, the locations of dynamic phase transition points change significantly 
depending upon amplitude and period of the external magnetic field as well as upon the active concentration of A- 
type components. Much effort has also been dedicated to clarify the hysteresis tools, such as coercivity, dynamic loop 
area as well as dynamic correlations between time dependent magnetizations and external time dependent applied 
field as a functions of period and amplitude of field as well as active concentration of of A- type components, and 
outstanding physical findings have been reported in order to better understand the dynamic process underlying present 
system. 
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1. Introduction 

It was shown, for the first time, by theoretically that when a ferromagnetic material with coupling J is exposed 
to a time dependent driving oscillating magnetic field, the system may not respond to the external magnetic field 
instantaneously, which gives rise to the existence of unusual and interesting behaviors due to the competing time 
scales of the relaxation behavior of the system and period of the external magnetic field (ll]. A typical ferromagnet 
exists in dynamically disordered (P) phase where the time dependent magnetization oscillates around value of zero 
for the high temperature and amplitude of field regimes. At this stage, the time dependent magnetization of system 
is capable to follow the external field with relatively small phase lag. However, it oscillates around a non-zero value 
which indicates a dynamically ordered (F) phase for low temperatures and small amplitude of field regimes. After 
that, a great deal of studies concerning the non-equilibrium phase transitions and also hysteresis behaviors of different 
types of magnetic systems has been investigated by using both experimental and theoretical techniques. For the sake 
of completeness, it is beneficial to talk about some of the prominent experimental and theoretical studies, respectively. 
Actually, most of the experimental works have been devoted to the dynamic phase transitions and hysteretic treatments 
of different types of thin films ESSiSIl. For example, in Ref. ||3l, dynamic magnetization reversal behavior 
of polycrystalline Ni 8 oFe 2 o films has been studied by applying a magnetic field along the easy axis of the sample 
and it is found that the hysteresis loop area A ^MdH^ is found to follow the scaling relation A oc HqQPT^'*' with 
a » 0.9, j3 0.8 and y - 0.38, where Hq, Q. and T are the amplitude, frequency of the external magnetic field and 
temperature, respectively. Moreover, it has been shown recently by Berger et al. in Ref. ||3l that uniaxial Co film 
sample under the influence of both bias (namely time independent) and time dependent oscillating magnetic fields in 
the vicinity of dynamic phase transition displays transient behavior for t < Tc, where r and are period and critical 
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period of the external applied field. Based on the experimental investigations we mentioned briefly above, it has been 
discovered that experimental non-equilibrium dynamics of considered real magnetic systems strongly resemble the 
dynamic behavior predicted from theoretical calculations of a kinetic Ising model. 

On the other hand, from the theoretical point of view, dynamic phase transitions as well as hysteresis properties 
of both bulk and finite size lattice systems have been clarified within the frameworks of Mean-Field Theory (MFT) 
iSS I3> 11.12.121, Effective-Field Theory (EFT) with single-site correlations 1 14,^. 
well as method of Monte Carlo (MC) simulations I 2 II 22 . 23 . 24, 25, 2^ 27, 28, 2S .Tic . 


m 20, M, as 
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331] . For example. 


by benefiting from a detailed large-scale MC simulation, in Ref. OOll . Park and Pleimling considered kinetic Ising 
models with surfaces subjected to a periodic oscillating magnetic field to probe the role of surfaces at dynamic phase 
transitions. They reported that the non-equilibrium surface universality class differs from that of the equilibrium 
system, although the same universality class prevails for the corresponding bulk systems. Moreover, by making use 
of MC simulation with local spin update Metropolis scheme, dynamic phase transition features and stationary-state 
behavior of a ferrimagnetic nanoparticle system with core-shell structure have been recently analyzed by some of us 
in Ref. IH. It has been observed that the particle may exhibit a phase transition from P to F phase with increasing 
ferromagnetic shell thickness in the presence of ultra-fast switching fields. Apart from these, there exists a limited 
number of theoretical studies about the non-equilibrium phase transition properties of magnetic systems containing 
quenched randomness resulting from random interactions between the spins with the same magnitudes or from a 
random dilution of the magnetic ions with non-magnetic species on the magnetic materials 118, Bl iilslIH. It is 


clear that these studies have a crucial role to have a better insight of the physics behind on real materials since a lot 
of magnetic materials have some small defects, and magnetic properties, i.e. phase transition temperature point of 
sample varies significantly depending on the type of defects. 

In this letter, for the first time we intend to determine the magnetic phase transition and hysteretic properties of 
a disordered binary ferromagnetic alloy of the type ApBi^p system under a magnetic field that oscillates in time. 
Such a quenched disordered system consists of two different species of magnetic components, namely A and B, and 
we select the magnetic components A and B to be as 5 a = 1/2 and S b = 1, respectively. Here, p refers to the 
concentration of type-A magnetic ion. The square lattice sites are randomly occupied by A and B atoms depending on 
the selected concentrations of magnetic components. Before going further, we should note that equilibrium or static 
features of such of disordered binary magnetic systems have been analyzed by means of several types of frameworks 
such as Perturbation Theor y ll^ MFT |37, 38, ^ 4^, Bethe-Peirls Approximati on jl^ . EFT with single-site 
coiTelations |42, 4^ £7, 4§] and method of MC simulation ||4§,l53, 51, 52, 531]. For instance, the critical 

properties of random mixtures of ferromagnetic and anti-ferromagnetic spin-spin interactions have been studied with 
MC simulation on a simple cubic lattice in Ref. It has been found that the system exhibits spin-glass phase 

characterized by a cusp-like peak in the susceptibility. Additionally, it is claimed by the author in Ref. 1391] that 
random-site binary ferromagnetic Ising model shows seven topologically different types of phase diagrams including 
a variety of multi-critical points within the framework MFT. 

The main motivation of the present study is to look answers for the physical facts underlying below questions: 


• What is the effect of the amplitude and frequency of the external magnetic field on the dynamic phase transition 
properties (i.e. critical temperature) of the considered system? 

• What kind of physical relationships exist between the magnetic properties (i.e. dynamic loop area, dynamic 
correlation and coercivity) of the system and the concentration of A-type magnetic components? In other 
words, whether the macroscopic treatments of the observable quantities of system depend on the concentration 
of A-type magnetic components or not. 

The plan of the remainder parts of the paper is as follows: In section (|2ll we briefly present our model. The results 
and discussions are presented in section (O, and finally section (|4]i contains our conclusions. 


2. Formulation 

The system simulated here, which is schematically shown in Fig. [1] is a disordered binary ferromagnetic alloy of 
the type ApBi^p which is defined on a two dimensional regular square lattice under the existence of a time dependent 
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driving oscillating magnetic field. The lattice sites are randomly occupied by two different species of magnetic com¬ 
ponents A and B with the concentration p and I - p, respectively. Thus, the Hamiltonian of the considered system can 
be written in the following form: 

+ (1 - d,)( 1 - 6j)S iS j +Sid- 6j)aiS ^ + (1 - 6i)6jS i(r] - //(O ^ ^ > 

(i.i) i 

where 7 > 0 is the ferromagnetic exchange interaction energy between i and j sites. The cr and S are conventional 
Ising spin variables which can take values of cr = +1/2 and S = +1,0 for the magnetic A- and B- components of 
system, respectively. 5, is a random variable which can take value of unity or zero depending on whether the site-/ 
is occupied by A or B ion, respectively. The first summation in Eq. ([T]l is over the nearest-neighbor pairs while the 
second one is over all lattice sites in the system. Hit) denotes the time dependent oscillating magnetic field described 
as H{t) - hosm(a>t), where t is time, ho and a> are amplitude and angular frequency of the driving magnetic field, 
respectively. The period of the oscillating magnetic field is given by t = ItiIlo. 



Figure 1: Schematic representation of the system of the type ApB{-p for five selected active concentration values of A-type magnetic component, 
namely (a) p = 0.0, (b) p = 0.25, (c) p = 0.5, (d) p = 0.75 and (e) p = 1.0. The blue and red colors correspond to A and B types of magnetic 
components which are randomly distributed. 


We simulate the system specified by the Hamiltonian in Eq. ([TJ on a L x L, where L — 128, square lattice 
under periodic boundary conditions applied in all directions by benefiting from MC simulation with single-spin flip 
Metropolis algorithm |3 3- The general simulation procedure we use in this study is as follows: The simulation 
begins from high temperature using random initial conditions, and then the system is slowly cooled down with the 
reduced temperature steps ksATIJ = 0.015, where configurations were generated by selecting the sites randomly 
through the lattice and making single-spin-flip attempts, which were accepted or rejected according to the Metropolis 
algorithm. The numerical data were generated over 50 independent sample realizations by running the simulations for 
20000 MC steps per site after discarding the first 10000 steps. This amount of transient steps is found to be sufficient 
for thermalization for the whole range of the parameter sets. 

Our program calculates the instantaneous values of the magnetizations ATiCf) and Msit), and also the total mag¬ 
netization Mrit) at time t as follows: 
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Mrit) = MA{t) + MbH), 


( 2 ) 


where Na denotes the total number of A ions (Na - pL^ while Nb represents the total number of B ions (Nb - 
(1 - p)LP')- Eurther, utilizing the instantaneous magnetizations stated in Eq. (|2l), we determine the dynamic loop area 
which measures the energy loss of the system due to the hysteresis, and dynamic correlations between time dependent 
magnetic field and magnetizations, respectively as follows 


Aa^-d) Ma{t)dH, 


C -1 

T 


Ma{t)Hit)dt, 


(3) 


here a - A,B and T. Namely, Aa,Ab, and Aj correspond to the dynamic loop areas while Ca, Cb, and Ct refer to the 
dynamic correlations of magnetic components A, B and total system, respectively. 
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In order to specify the location of dynamic phase transition point at which F and P phases separate from each 
other, we use and check the thermal variation of dynamic heat capacity and which is defined as; 


C = 


dEt 


Coop 


dT 


(4) 


where T is temperature while Ecoop is the cooperative part of the energy over a full cycle of the external applied 
magnetic field which is defined according to the following equation: 


-‘Coop 


1 

tL?- 


jJ] + (1 - 5/)(l - ^j)SiSj + Sid - 6j)(TiSj + (1 - dddjSitTj] 


(Uj) 


dt. 


(5) 


3. Results and Discussion 

In this section, in order to explain the dynamic evolution of the magnetic system in detail, we will focus our atten¬ 
tion on non-equilibrium phase transition properties as well as stationary-state treatment of the quenched disordered 
binary ferromagnetic alloy system under a time dependent oscillating magnetic field. We argue and discuss how the 
amplitude and period of the external oscillating magnetic field as well as concentration of A-type magnetic component 
affect the dynamic critical nature of the system. We will also shed light on the physical facts lying behind the thermal 
dependencies of dynamic loop areas as well as correlations between time dependent magnetizations and magnetic 
field for some selected combination of system parameters. As a final investigation, we examine the variations of 
coercivity and also hysteresis behavior as functions of system parameters, as well as active concentration of A-type 
magnetic component. 

Before going further, it should be emphasized that for HqIJ- 0.0, i.e. static case, our numeric MC simulation 
findings are completely in accordance with the recently published work 0531] where thermal equilibrium phase transi¬ 
tion properties of the system we consider have been addressed by utilizing MC simulation with Metropolis algorithm. 
Let us start to discuss the non-equilibrium phase transition properties of the system. In Figs. |2la-c), we give the global 
dynamic phase boundaries that distinguish between F and P phases in a (p - ksTclJ) space with five reduced external 
field amplitude values, such as ho/J - 0.0,0.25,0.5,0.75 and 1.0 for some considered values of oscillation periods, 
namely r = 50,100 and 200. At first sight, one can easily see from the figure that, (i) for a fixed set of values of ho/J 
and T, when the concentration of spin-1/2 atom is increased starting from p - 0.0 11561] . F to P phase transition point 
moves to a lower value in temperature axis. This is because of the fact that the system tends to become disordered 
due to the occurrence of two different types of randomly occupied magnetic components as the active concentration 
of spin-1/2 is increased starting from p - 0.0. 

It is possible to say that as the p value is increased then the energy contribution which comes from spin-spin 
interactions gets smaller. Consequently, the F to P phase transition point moves to a lower temperature, due to the en¬ 
ergy resulting from the temperature and (or) magnetic field which overcome the ferromagnetic spin-spin interactions. 
Therefore, the F regions in (p - ksTc/J) plane get narrower. Moreover, {if) for a fixed set of values of p and r, we 
see that when ho/J increases, then the magnetic energy supplied by the external applied field dominates against the 
ferromagnetic spin-spin exchange interactions, and hereby, the studied system can relax within the oscillation period 
T of the external field which gives rise to a reduction in the transition temperature. In order to see this physical fact, 
one can compare any two values of applied field amplitudes, for example ho/J - 0.25 with 0.5 for p = 0.8 in Figl2ja). 
Another important finding of our numerical simulation is that the location of the phase transition temperature sensi¬ 
tively depends on the applied field period. (Hi) For a fixed set of values p and ho/J, it is obvious that as r is increased, 
the phase transition point is lowered since decreasing field frequency gives rise to a decreasing phase delay between 
the magnetizations and magnetic field (i.e. the magnetizations can follow the oscillating driving magnetic field) and 
this makes the occurrence of the dynamic phase transition easy. The physical discussions mentioned above can be 
easily seen by checking any two values of oscillation periods such as t = 100 and 200 for fixed values of ho/J - 0.5 
and p = 0.8. 

We draw in Figs. |3ja-c) the active concentration of A-type magnetic component, reduced amplitude and oscillation 
period of field dependencies of dynamic correlation versus temperature curves. One can able to deduce from the fig¬ 
ures that, at relatively low temperature regions where the system exists in F phase, the dynamic correlations between 
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Figure 2: Dynamic phase diagrams in (p — kgTclJ) planes obtained from the peaks of the dynamic heat capacity for values of Hq/J = 
0.0,0.25,0.5,0.75 and 1.0. The curves are demonstrated for three values of oscillation periods: (a) t = 50, (b) r = 100 and (c) r = 200. 
The characters F and P refer to the dynamically ferromagnetic and paramagnetic phases, respectively. 
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Figure 3: Thermal variations of dynamic total correlations corresponding to the dynamic phase diagrams shown in Figs.|2la-c). The parameters we 
study are remarked on the graphs. 


time dependent magnetizations of the system and forcing field are close to zero since thermal energy is almost negli¬ 
gible, and the ferromagnetic spin-spin interactions are dominant against the field energy. Moreover, the system begins 
to show a shallow dip behavior, which may become negative depending on the system parameters, in the vicinity of 
dynamic phase transition point of system. For example, the dynamic correlation curves exhibit concentration, reduced 
amplitude and oscillating period of held induced negative shallow dip behavior in Figs. [3 a), (b) and (c), respectively. 

In the following analysis, in order to show what happens in dynamic loop area versus temperature curves at 
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Figure 4: Temperature variations of dynamic total loop areas for varying system parameters which are indicated in related graphs. 


varying system parameters, we give in Figs. SJa-c) thermal variations of dynamic loop areas corresponding to the 
dynamic phase diagrams illustrated in Figs. 12a-c). As indicated in section (|2]l, dynamic loop area characterizes 
energy dissipation due to the hysteresis. It can be readily seen from the figures that the dynamic loop area curves 
obtained for a wide range of system parameters reveal a smooth and rounded cusp located at a certain temperature 
above the dynamic phase transition temperature of system. It would be interesting to see the origin of the physical 
mechanisms underlying on the thermal treatment of dynamic loop areas. Keeping this in mind, we separated the 
temperature space into five stages which are shown in Fig. |5]where dynamic loop areas of two different parts as well 
as of overall system have been plotted for selected values of system parameters such as p = 0.5, hojj - 0.5 and 
T = 100. At stage-I, the temperature is sufficiently low so that both magnetic components A and B of the system 
show a dynamic ferromagnetism. In other words, the time dependent magnetizations of the system can not follow 
the external field simultaneously. Hence, dynamic loop areas of the system are close to zero. If we look at stage-II, 
where the location of the stage-II corresponds to the dynamic critical temperature of the system for the parameters we 
select, we can see that an increment in temperature gives rise to a reduction of ferromagnetism, and time dependent 
magnetizations of the system begin to relatively respond to the time dependent external magnetic field. As a result of 
this, hysteresis loop areas of both magnetic components A and B of the system get wider. Before going further, we 
note that even though the stage-II is a phase transition point, the dynamic loop areas does not show a maximum peak 
at the position of stage-II. This is because of such MC computer experiments include of thermal fluctuations of spins. 
If one considers the same problem with MFT (or EFT) which neglects (or partially includes) the thermal fluctuations 
of spins, and focuses on the thermal evolution of dynamic loop areas of the system, with a high probability, one 
can see that thermal variation of dynamic loop area may exhibit a maximum at the dynamic phase transition point 
ifl^l^ . Let us continue with stage-III where the dynamic loop areas of the system present a maximum peak, both 
magnetic components A and B of the system show a dynamic paramagnetism, i.e., the time dependent magnetizations 
of the system can follow the external field with a relatively small phase lag. As the temperature is increased above 
the dynamic transition point, we reach the stages-IV and -V, where the value of dynamic loop areas begin to reduce, 
respectively. Apart from the changing values of the dynamic loop area corresponding to different stages we determine, 
shapes of hysteresis curves sensitively depend on the selected system parameters which will be discussed below. 

In order to understand the varying temperature influences on the hysteresis treatments of the system, we give 
in Fig. |6]the hysteresis curves for both magnetic components A, B and overall of the system for selected system 
parameters corresponding to the diagram depicted in Fig. |5] When we look at the hysteresis curves obtained for 
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Figure 5: Thermal variations of dynamic loop areas of magnetic components A, B and overall system for considered values of system parameters, 
such as p = 0.5, h^jj = 0.5 and r = 100. The labeled stage-i {i =1, II,..,V) refers to the five different temperature locations indicated in graph. 


stage-I, we see that hysteresis loops are asymmetric around zero value indicating an F phase. If one moves throughout 
stage-I —> stage-II way, it can be easily seen that the asymmetric shapes of hysteresis curves begin to slowly break 
down owing to the existence of a higher thermal energy than before. When the studied system attains the stage-III, 
it is possible to observe an example of square-like shape of hysteresis curves. In this zone, the system completely 
exhibits a P character. As the temperature increased, we reach the stage-IV, and we see that shape of hysteresis loops 
begin to change from square-like to sigmoidal-like shaped due to the thermal agitations. In addition to these, further 
increment in temperature leads to the occurrence of a narrower sigmoidal-like hysteresis curves which are explicitly 
shown in stage-V of figure. Similar of observations have been recently found in ultrathin Blume-Capel films under a 
magnetic field which oscillates in time 15711 . 

In the following two analysis, namely in Figs. |2]and|8] we investigate the oscillation period, reduced amplitude 
of external applied field and concentration of A-type magnetic component dependencies of the coercivity as well 
as hysteresis curves of the system. The aforementioned physical quantities have been calculated at a temperature 
T - 0.87’°(p), where T^(p) is the static phase transition point in the absence of the external field, and it depends on 
the studied active concentrations of magnetic atoms. The reason why we choose such a temperature is that the system 
can be capable of undergoing a purely mechanical phase transition, such as applied field period or reduced amplitude 
induced phase transition. The numerical data were collected for 100 cycles of the external field after discarding the 
first 100 cycles of field have been discarded to obtain a stationary state behavior. 

We give the period of external field dependencies of the coercivity of the system in Fig. 0 with varying values 
of reduced applied field amplitudes, such as ho/J - 0.25,0.5,1.0,2.0 and 4.0 with three considered values of active 
concentration of A-type magnetic component, namely p - 0.0, 0.5 and 1.0. It can be said by focusing on Fig. [Tja) that 
depending on the reduced applied field amplitude, it is possible to observe a large coercive field at low applied field 
periods. Moreover, it is found that coercivity curves of system exhibit a sudden variation with increasing applied field 
period whereas for sufficiently high t they exhibit a stable profile. Another important point we want to emphasize is 
that the system has no coercivity value in the case of low values of applied field period and reduced amplitude, for 
example t = 100 and Iiq/J = 0.25. In this regime, the system exists in F phase. However, with a significant increment 
in value of amplitude for the same period of field causes the existence of a purely mechanical reduced amplitude 
induced phase transition, and coercivity treatment begins to show itself naturally. The nearly same discussions are 
also valid for Figs|7]^b-c) which are plotted for p - 0.5 and p - 1.0, respectively. As the type-A magnetic component 
is introduced into the lattice, the coercivity value decreases prominently since the energy originating from the time 
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Figure 6: Hysteresis curve evolutions of both magnetic components A and B and also overall of system at different five temperature locations for 
the same system parameters with figure [5] 


dependent forcing field overcomes the ferromagnetic spin-spin interaction term between spin pairs in the lattice. 
Moreover, we should note that at the relatively low values of the amplitude and the period of the external field, even 
though some small number of spin flip occurs, the considered alloy system shows a dynamic ferromagnetic character. 
The physics mentioned briefly here is valid for all concentration values of the system. 
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Figure 7: Variations of the coercive fields as a function of the oscillation period for several values of ho/J with (a) p = 0.0, (b) p = 0.5 and (c) 
p = 1.0, respectively. 


As a final investigation, we present the evolution of some typical hysteresis curves corresponding to the coercivity 
properties illustrated in Fig. [T] Throughout Figs. m a-c), we fix the active concentration of magnetic components 
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Figure 8: Influences of (a) reduced amplitudes for r = 100, (b) relatively low periods, and (c) high periods of field for ho/J = 0.5 on the hysteresis 
loops of the system corresponding to the figure|2jb). 


as /? = 0.5 and change the external field components, namely period and reduced amplitude of field. In Fig. [Ha), 
we select the period of field as t = 100 to investigate the influences of varying reduced amplitudes of field on 
hysteresis curves. It is clear from the figure that while the system exists in dynamically ordered phase and hysteresis 
curve is asymmetric around zero value at the relatively low reduced external field amplitude, i.e. holJ- 0.25, with 
increasing amount of ho/J the system starts to show a dynamically paramagnetic character, and hysteresis curve is 
symmetric around zero value. According to our numerical findings, another important observation is that the shape of 
the hysteresis curves sensitively depend on the applied field amplitude, for example, while the hysteresis loop is broad 
square-like shaped for value of ha/J - 0.5, it is sigmoidal-like shaped for values of ho/J — 1.0, 2.0 and 4.0. In Fig. 
[HJb), we give the dynamic hysteresis loops obtained for selected values of applied field period such as t = 50,100,200 
and 400 with fixed value of ho/J = 0.5. At first glance, by comparing the curves obtained for t = 50 and 100, one 
can clearly deduce from the figures that (i) remanent magnetization and also coercivity values of the system tend to 
change, (ii) the hysteresis curves begin to widen indicating a period of field induced phase transition with increasing 
T. If we increase t, the hysteresis curves becomes broad square-like shaped. On the other hand, in Fig. Htc), it is 
possible to see that the hysteresis curves are sigmoidal shaped, and the width of the loops becomes narrower but does 
not vanish at the relatively high values of applied field period regimes. This type of behavior has been reported in 
M where non-equilibrium phase transition properties of a nanocube system is investigated with MC method of 
simulation. 


4. Conclusion 

In conclusion, by making use of MC computer experiment method with Metropolis algorithm, we have looked 
for answers for the questions discussed in section [T] As we mentioned above, as a disordered binary ferromagnetic 
alloy with coupling J, we have selected a system consisting of spin-1/2 and spin-1 components, where the spin-1/2 
and spin-1 components are distributed randomly depending upon the controllable concentration value. The situations 
of p = 1.0 and 0.0 correspond to the completely pure kinetic Ising and Blume-Capel models without single-ion 
anisotropy, respectively. Namely, for both two values of p we remark, there is no lattice including disorder originating 
from different species of components in the system. However, it is obvious that except from these two values of p, 
the magnetic system composes of randomly distributed spin-1/2 and spin-1 components where the weights of the 
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components strongly depend on the concentration values. After detailed numerical operations, the most prominent 
findings underlined in the present paper can be briefly summarized as follows; 

• For a fixed set of values of h^jj and t, we found an almost linear decrease the critical temperature with increas¬ 
ing type-A content. 

• We see that when holJ increases, then the magnetic energy supplied by the external applied field dominates 
against the ferromagnetic spin-spin exchange interactions, and hereby, the studied system can relax within the 
oscillation period t of the external field which gives rise to a reduction in the transition temperature, for a fixed 
set of values p and t. 

• For a fixed set of values p and hojj, as r is increased, the phase transition point is lowered because of increasing 
field frequency gives rise to a growing phase delay between the time dependent magnetizations and forcing 
applied field. 

The situations corresponding to the explanations we emphasize above can be easily seen by looking the dynamic phase 
boundaries separating the F and P phases which have been constructed in active concentration of spin-1/2 atoms and 
dynamic critical temperature planes at various Hamiltonian parameters in Fig. [T] Our simulation results also suggest 
that: 

• For studied values of system parameters, thermal variations of dynamic correlations exhibit period and reduced 
amplitude of field as well as active concentration of A-type magnetic component induced local dip behavior. 

• The locations of maximum lossy points which are characterized by hysteresis loop area show a decreasing 
tendency in temperature value with increasing p values starting from zero for considered values of Hamiltonian 
parameters. Other system parameter dependencies of dynamic loop area has also been discussed in the present 
study. 

• Applied field period dependencies of the coercivities obtained for three concentration values of p with varying 
reduced amplitudes of field reveal that they exhibit a stable profile at the relatively high period of field whereas 
they show a sudden change in the low period of field. It is also found that for a fixed set of ho jj and t, as the p 
value increases starting from zero, the corresponding coercivity value begins to decrease this is because of the 
amount of spin-1/2 atoms increases in system. 

• Shape of the hysteresis curve of the system strongly depends on the taking Hamiltonian parameters into consid¬ 
eration. For example, while it is possible to observe an example of a broad square-like shaped hysteresis curve 
at the relatively low period of field, the shape of it begins to change to be a sigmoidal-like shape with increasing 
period of field. 


All of the observations found in this work show that the amplitude and frequency of the driving time dependent 
magnetic field as well as active concentration of A-type magnetic component have an important influence on the 
thermal and magnetic properties, such as coercivity and critical temperature of the system. 

Finally, we note that such a binary alloy system can actually contain four different values of relative A-B interac¬ 
tion energy, namely Jaa, Jbb, Jab, Jba corresponding to the randomly distributed components between A - A, B - B, 
A-B and B - A, respectively. According to the findings of the previously published equilibrium phase transition 
studies mentioned in section [T] disordered binary alloy system with different signs and unequal magnitudes of the 
spin-spin interactions and single ion-anisotropy reveals interesting thermal and magnetic behaviors such as existence 
of re-entrant, compensation and also spin-glass behaviors. In this work, we have focused on only the non-equilibrium 
phase transition properties of disordered binary alloy system in cases of Jaa - Jbb - Jab - Jba - J and also zero 
single-ion anisotropy. So, determination of the non-equilibrium phase transitions features of a disordered binary alloy 
including different types of spin-spin interactions and also single-ion anisotropy may be an interesting study. Fur¬ 
thermore, based on the earlier systems 1 1^ 19, 35], it is possible to categorize the frequency dispersions of dynamic 
loop areas into three distinct types depending on the system parameters. Keeping in this mind, it would also be very 
interesting to see the frequency dispersions of dynamic loop areas of the same system. 
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